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Abstract 

We classify integrable third order equations in 2 + 1 dimensions which generalize the 
examples of Kadomtsev-Petviashvili, Veselov-Novikov and Harry Dym equations. Our ap- 
proach is based on the observation that dispersionless limits of integrable systems in 2 + 1 
dimensions possess infinitely many multi-phase solutions coming from the so-called hydrody- 
namic reductions. In this paper we adopt a novel perturbative approach to the classification 
problem. Based on the method of hydrodynamic reductions, we first classify integrable quasi- 
linear systems which may (potentially) occur as dispersionless limits of soliton equations in 
2 + 1 dimensions. To reconstruct dispersive deformations, we require that all hydrodynamic 
reductions of the dispersionless limit arc inherited by the corresponding dispersive counter- 
part. This procedure leads to a complete list of integrable third order equations, some of 
which are apparently new. 

MSC: 35L40, 35Q51, 35Q58, 37K10, 37K55. 
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1 Introduction 



The classification of integrable systems has been a topic of active research from the very begin- 
ning of soliton theory. In 1 + 1 dimensions, this resulted in extensive lists of integrable equations 
within particularly important subclasses [23J, which were obtained by means of the symmetry 
approach. Although this technique generalizes to 2 + 1 dimensions, one encounters additional 
difficulties due to the appearance of non-local variables [21] . A way to bypass the problem of 
non- locality, known as the perturbative symmetry approach |25j . provides an efficient way to 
classify soliton equations in 2 + 1 dimensions. In this framework, one starts with a linear equa- 
tion having degenerate dispersion law [39], and reconstructs the allowed nonlinearity. However, 
few classification results have been obtained so far. In fact, most of the (2 + l)-dimensional 
examples known to date were derived by postulating a special structure of the corresponding 
Lax pair, see e.g. [50], [55] . 

In this paper we adopt a novel approach to the problem of classification of scalar third order 
soliton equations in 2 + 1 dimensions with the 'simplest' possible non-localities, 

ut = F(u, w, Du, Dw), 

and 

u t = F(u, v, w, Du, Dv, Dw), 

respectively. Here u(x,y,t) is a scalar field, and the non-local variables v(x,y,t) and w(x,y,t) 
are defined via w x = u y and v y = u x , equivalently, w = D~ D y u, v = D~ l D x u. The symbols 
Du, Dv, Dw denote the collection of all partial derivatives of u, v, w with respect to x and y up 
to the third order. In fact, it is sufficient to allow only y-derivatives of w and x-derivatives of v. 
We will refer to the above equations as the 'non-symmetric' and 'symmetric' cases, respectively. 
We assume that in both cases the dependence of the right hand side F on the derivatives of u 
and w (resp, u, v, w) is polynomial, where the coefficients are allowed to be arbitrary functions 
of u and w (resp, u,v,w). Explicitly, in the non-symmetric case we have 

u t = (pu x + ipu y + r]w y + e(...) + e 2 (...), w x = u y , (1) 

where ip,ip,r] are functions of u and w, while the terms at e and e 2 are assumed to be homogeneous 
differential polynomials of the order two and three in the derivatives of u and w, whose coefficients 
can be arbitrary functions of u and w. We use the following weighting scheme: u and w are 
assumed to have order zero, their derivatives u x ,u y ,w x ,w y are of order one, the expressions 
u xx , u xy , u yy ,w yy ,u x ,u x u y , u y ,u x w y , u y w y , w 2 are of order two, etc. Thus, the term at e is a linear 
combination of the ten second order expressions whose coefficients can be arbitrary functions of 
u and w. The most familiar example within the class ([T]) is the Kadomtsev-Petviashvili (KP) 
equation, 

u t = uu x + w y + e 2 u xxx , w x = u y . 
Similarly, in the symmetric case we consider equations of the form 

ut = <pu x + ipiiy + rjWy + tv x + e(...) + e 2 (...), w x = u y , v y = u x , (2) 

here (p,tp,rj,T are functions of u, v and w. A canonical example of the form ([2]) is the Veselov- 
Novikov (VN) equation, 

U t = (UV) X + (Uw)y + C 2 (U XXX + Uyyy), W X = Uy, Vy = U X . 
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In Sect. 2 we bring together other known examples of the form ([T]) and which include the 
KP, VN, Harry Dym equations and their modifications. 

Our approach to the classification problem is based on the following key observations. 

• Dispersionless limits of integrable soliton equations in 2 + 1 dimensions possess 
infinitely many hydrodynamic reductions. 

In particular, dispersionless limits of Eqs. (fl|) and (J2]), 

U t = (PU X + l{jUy + TJWy, W X = Uy, (3) 

and 

U t = (PU X + IpUy + TjWy + TV X , W X = Uy, Vy = U X , (4) 

should possess infinitely many hydrodynamic reductions and, thus, must be integrable in the 
sense of [10]. It was observed in [10] that the method of hydrodynamic reductions provides 
an efficient classification criterion. Thus, as a first step, in Sect. 3 we classify integrable first 
order equations of the form ([3|) and which may (potentially) occur as dispersionless limits 
of integrable equations of the form ([T]) and ([2]). We emphasize that the requirement of being a 
dispersionless limit of a third order soliton equation imposes further severe constraints, so that 
very few particular cases obtained in Sect. 3 do actually survive. 

Given an integrable dispersionless equation, one needs to reconstruct dispersive deformations. 
In 1 + 1 dimensions, this problem has been a subject of extensive research in [3 El El E] , see also 
[lj. In 2 + 1 dimensions, the reconstruction procedure is based on the following key observation 

m- 

• Hydrodynamic reductions of dispersionless limits of integrable soliton equa- 
tions can be deformed into reductions of the corresponding dispersive coun- 
terparts (strictly speaking, this is only true if the dispersionless limit is linearly 
non-degenerate, see Sect. 4). Furthermore, the requirement of the inheritance 
of all hydrodynamic reductions allows one to efficiently reconstruct dispersive 
terms in 2 + 1 dimensions. 

This suggests the following alternative definition of the integrability: 

A (2 + 1) -dimensional system is said to be integrable if all hydrodynamic reductions of its disper- 
sionless limit (which is assumed to be linearly non-degenerate) can be deformed into reductions 
of the corresponding dispersive counterpart. 

Although this property is satisfied for all known integrable equations whose dispersionless limit 
is not totally linearly degenerate, it would be important to formulate more precise statements 
about the equivalence of our definition with more 'conventional' approaches to the integrability. 

The procedure of the reconstruction of dispersive terms is thoroughly illustrated in Sect. 4, 
where we examine case- by-case all integrable dispersionless limits from Sect. 3. Our calculations 
result in a complete list of integrable (2+l)-dimensional equations, some of which are apparently 
new. It is important to emphasize that, although our approach is based on the requirement of 
the inheritance of hydrodynamic reductions, all examples from the final list do actually possess 
conventional Lax pairs. Altogether, we found three new equations. One of them is 

Ut = (/3w + (3 2 u 2 )u x — 3(3uu y + w y + e 2 [B 3 (u) — j3u x B 2 {u)}, w x = u y , (5) 
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here B = (3uD x — D y , /3=const. It possesses the Lax pair 

i> xy = (3uil) xx + ^ V', 

1p t = e 2 (3 3 U 3r lp xxx - e 2 4>yyy + 3e 2 P 2 UUyl/j xx + f3wip x . 

The second example is 

ut = -f5 2 u 3 u x + (w — 3(3u 2 )u y + uw y + e 2 [B 3 (u) — f3u x B 2 (u)}, w x = u y , (6) 
3 

here again B = f3uD x — D y , /3=const. The corresponding Lax pair is 

1 



tpxy = f3uil) xx + ^Ul/l, 

f3 2 

tpt = e (3 3 u 3 i(} xxx - e 2 ip yyy + 3e 2 p 2 uu y ^ xx + — u 3 %(j x + wifjy + f3uuyip. 



We point out that similar Lax operators appeared in the context of the (2 + l)-dimensional 
Camassa-Holm equation [JO]. Our last example is a deformation of the Harry Dym (HD) equa- 
tion, 

U t = ^U X - 2WUy + UWy - — ( - J , WX =Uy, (7) 



u u \u, xxx 



for 5 = it reduces to the standard HD equation (Example 6 of Sect. 2.1). It has the Lax pair 
L t = [A, L] where 



L "^ 2 + 7!^ + 4^' 

All three examples belong to the non-symmetric case. In the symmetric case we have no new 
equations apart from those listed in Sect. 2.2. This leads to the following main result: 

Theorem 1 Equations ^ - ^ along with the known examples of KP, non- symmetric VN, 
HD equations and their modifications provide a complete list of integrable equations of the form 
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(OP with n ^ whose dispersionless limit is linearly nondegenerate: 

KP equation ut = uu x + w y +e 2 u xxx , 

mKP equation u t = (w — u 2 /2)u x + w y +e 2 u xxx , 

B 2 

Gardner equation ut = (8w — — u 2 + 8u)u x + w y +e 2 u xxx , 



VN equation ut = (uw)y+e 2 u yy y, 



. . 2 / 3 ul 
mVN equation u t = (uw) y +e I u yy — — — 

e 2 (\ 

HD equation u t = —2wu v + uw v — 

y y u \ u 



xxx 
2 



5 e z (\ 

deformed HD equation u t = —u x — 2wu y + uw y 



u° u \u 

2„.2 



Equation © iH = -P 2 u^u x + (w — ?>Bu 2 )u y + uw y + e 2 [B 3 (u) — 8u x B 2 (u)]. 



Equation © u t = (3w + 3 u )u x — 33uu y + w y + e [B (u) — 3u x B (u)], 

4 

In the symmetric case there exist only two examples of integrable equations of the form (0j with 
n,r^0: 



VN equation ut = (uv) x + (uw) y + e 2 u xxx + e 2 u 

( 3u 2 \ ( 
mVN equation ut = (uv) x + (uw) v + e 2 I u xx ) + e 2 

V 4 Wx V 



Zu 2 y 



4 u 



The proof is summarised in Sect. 4. Under the substitution w = 0, u y = the equations 
dH) reduce to 

u t = e 2 p 3 (u 3 u xxx + 3u 2 u x u xx ) + 3 2 u 2 u x 

and 

n t = e 2 3 3 (u 3 u xxx + 3u 2 u x u xx ) + -^n 3 ^, 

respectively. In this form, they have appeared in [M], see also [35] and references therein. It 
was pointed out (see e.g. [HHIEHIEE]) that there exist differential substitutions bringing these 
equations to a constant separant form. It would be interesting to find out whether Eqs. (|5]) - 
([7|) are related to any of the known soliton hierarchies: the main problem here is that the above 
differential substitutions do not extend to 2 + 1 dimensions in any obvious way. 
Remark 1. The examples of VN and mVN equations show that different (2 + l)-dimensional 
equations may have one and the same dispersionless limit. 

Remark 2. Our approach to the classification problem does not apply to non-symmetric 
equations with n = (or symmetric equations with n = r = 0). As we explain in Sect. 3, 
these conditions are equivalent to the reducibility of the dispersion relations of the corresponding 
systems ©, @. A familiar example within this class is the so-called 'breaking soliton' equation, 

u t = 2wu x + 4uu y - e 2 u xxy , w x = u y , 

see e.g. [3]. Here tp = 2w, tp = 4u, n = 0. Equations of this type are not amenable to the 
method of hydrodynamic reductions, and require an alternative approach. 
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2 Known Examples 

2.1 Non-symmetric case 

Here we bring together known examples of soliton equations whose dispersionless limit is of the 
form ([3]). The relation w x = u y will be automatically assumed whenever w appears explicitly in 
the equation. Examples 1-6 list third order equations. Examples 7-10 correspond to equations 
of order five, or differential-difference equations. 

Example 1. The Kadomtsev-Petviashvili (KP) equation, 

ut = uu x + w y + e 2 u XXXl (8) 

arises in mathematical physics as a two-dimensional generalization of the KdV equation. Its 
dispersionless limit (dKP equation), 

u t = uu x + w y , (9) 

also known as the Khokhlov-Zabolotskaya equation [37], is of interest in its own, playing impor- 
tant role in non-linear acoustics, gas dynamics and differential geometry. 

Example 2. The modified KP (mKP) equation, 

u t = (w - u 2 /2)u x + Wy + e 2 u xxx , (10) 

has the dispersionless limit 

ut = (w - u 2 /2)u x + w y . (11) 

Example 3. The (2 + l)-dimensional version of the Gardner equation is of the form [20J, 

(3 2 

u t = (flw - —u 2 + 5u)u x + w y + e 2 u xxx , (12) 



which reduces to the KP or mKP equations upon setting (3 = or 5 = 0, respectively. Its 
dispersionless limit has the form 

(3 2 

u t = (j3w — —u 2 + 5u)u x + w y . (13) 
Example 4. The non-symmetric version of the Veselov-Novikov equation [32 1 \27 \ 13] . 

U t = (Uw)y + e 2 Uyyy, (14) 

has the dispersionless limit 

u t = {uw) y . (15) 
Example 5. The non-symmetric version of the modified Veselov-Novikov equation [2], 

Ut = (UW)y +e 2 [ Uyy - - j , (16) 



u 
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has the same dispersionless limit as in the previous example, 

U t = (Uw)y. 



e 2 (1 



Example 6. The Harry Dym equation [20J, 

ut = —2wu y + UWy 
(set u = 1/u to obtain the equation from |20j), has the dispersionless limit 

Ut = —2WUy + UWy. 



Example 7. The fifth order version of the Harry Dym equation is 

e 2 e 4 / 1 

U t = SWUy + UWy l{u 2 U XX y ~ 3u(U X Uy) X + QU^Uy) + — I — 

see |20j . Its dispersionless limit has the form 

Ut = 2>WUy + UWy. 

Example 8. The Toda lattice is a system of two differential-difference equations 

eu t = u (w(y) - w(y - e)), 
ew x = u(y + e) - u(y), 

or 

e e 2 e n 

U t /U = Wy~ -Wyy + —Wyyy ( ~ 1 ) " + 1 — W n y + . . . , 
W X — Uy + -Uyy + —Uyyy + • • • + ~{ U Uy + 

Its dispersionless limit is 

Ut = UWy. 

Example 9. The nonlocal Toda lattice equation is 

cr(x-\-e,y-\-e) — a <j — <r(x — e,y — e) 

ea xt = e - e * , 

see [3T]. Its dispersionless limit is 

or, setting o x = u, a y = w, 

Ut = e U+W (U X + 2Uy+Wy). 

Example 10. The BKP and CKP equations are of the form 



i 



u t - 5(1* + w)u x - 5u Wx + 5w y + e (uu xxx + w xxx + u xxx ) - - u xxxxx = 0, 



and 

.2 , „..\„. r„.„.. , c... I ,2 



5 e 4 

«i - 5(n 2 + w)u x - 5uw x + bw y + e l {uu xxx + w xxx + -u xxx ) - — u xxxxx = 0, 
respectively [20J. Their dispersionless limits coincide: 

"y — UWy. 



u t = 5(n 2 + w)u x + buuy — 5w y . 



2.2 Symmetric case 

Here we list known examples of the form ([2]). The relations v y = u x and w x = u y will be 
automatically assumed whenever v and w appear explicitly in the equation. It is quite remarkable 
that the 'symmetric' list is very restrictive, and contains only two examples. 

Example 1. The Veselov-Novikov equation, 

ut = (uv) x + (uw) y + e 2 u xxx + e 2 Uyyy, (30) 

was introduced in [32], [27]. It has the dispersionless limit 

u t = (uv) x + (uw) y . (31) 



Example 2. The modified Veselov-Novikov equation, 

u t = {uv) x + {uw)y + e \^u xx - -— j + e \u yy - j , (32) 

was first introduced in [2] (in a somewhat different form). It has the same dispersionless limit 
as in the previous example, 

u t = (uv) x + (uw) y . (33) 



3 Classification of integrable dispersionless limits 

In this section we classify integrable dispersionless equations of the form @ and @j which 
may potentially occur as dispersionless limits of integrable soliton equations of the form ([T|) and 
(0), respectively. The integrability conditions are derived based on the method of hydrodynamic 
reductions. For the convenience of the reader, we briefly recall the main steps of this construction. 
As proposed in [10] , the method of hydrodynamic reductions applies to quasilinear equations of 
the following general form: 



A(u)u t + B{u)u x + C{u)u y = 0; (34) 

here u = (u , u 171 ) 1 is an m-component column vector of the dependent variables, and A, B, C 
are m x m matrices. The method of hydrodynamic reductions consists of seeking multi-phase 
solutions in the form 

u = u(R\...,R N ) (35) 

where the 'phases' R l (x,y,t) are required to satisfy a pair of consistent equations of hydrody- 
namic type, 

R y = fi l (R)R x , R\ = \ l (R)R x . 

We recall that the consistency conditions, R l yt = R\ y , imply the following restrictions for the 
characteristic speeds [i % and A 1 : 

djtf djX 1 
/m> - fi { ~ \i - A* ' 
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i / j, di = d/d R i, see [3D]. The substitution of the ansatz ([35]) into (f34"|) leads to a complicated 
over-determined system of PDEs for the functions u(R), ^ 1 {R) and \ % (R) whose coefficients 
depend on the matrix elements of A, B, C, and their derivatives. In particular, the characteristic 
speeds fi l (R) and X l (R) satisfy an algebraic relation det(XA + B + fj,C) = which is nothing but 
the dispersion relation of the system (|34p . We will assume that the dispersion relation defines 
an irreducible algebraic curve of degree m. 

Definition [TDJ- System (fff^| ) is said to be integrable if, for any number of phases N, it possesses 
infinitely many N -phase solutions parametrized by 2N arbitrary functions of one variable. 

The requirement of the existence of such solutions imposes strong constraints on the matrices 
A, B, C, which can be effectively computed. Although these constraints are quite formidable in 
general, there exists a simple necessary condition for the integrability which can be expressed 
in an invariant differential geometric form as follows. Let us first introduce the m x m matrix 

V = (aA + pB + jC)- 1 ^ + f3B + 7C) 

where a, f3, 7 and a, $, 7 are arbitrary constants. Given a (1, l)-tensor V = [u*-] 3 let us introduce 
the following objects: 
Nijenhuis tensor 

M) k = v]d u vv\ - v\d u vv) - v l p {d uJ v p k - d uk v?), 

Haantjes tensor 

One has the following result. 

Theorem 2 [12] The vanishing of the Haantjes tensor is a necessary condition for the integra- 
bility of the system (5$) . 

Since the Haantjes tensor can be obtained using computer algebra, one gets an efficient inte- 
grability test (notice that all components of the Haantjes tensor have to vanish for any values 
of the constants a,/3,7 and 7). These necessary conditions are very strong indeed, and 
in many cases turn out to be sufficient. We point out that, for m = 2, the Haantjes tensor 
vanishes identically and does not produce any non-trivial integrability conditions. In this case 
one proceeds as follows: let us multiply (I34|) by A^ 1 , and diagonalize B (this is always possible 
in the 2-component case). Thus, without any loss of generality one can assume 

In this particular normalization, the integrability conditions for 2x2 systems were obtained 
in [TT]. These conditions constitute a system of second order constraints for the coefficients 
a,b,p,q,r, s which can be easily tested. Let us now apply this approach to the classification of 
integrable systems of the form ([3j) and (jlj). 

3.1 Non-symmetric dispersionless limits 

Given an equation of the form (|3|), 

u t = (pu x + ipu y + rjWy, 
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let us first rewrite it in matrix form (|34|) as follows: 



This system is now in the form as studied in [TT]. The integrability conditions reduce to a 
system of second order partial differential equations for the coefficients (p, ip and rj, which can 
be derived from the general integrability conditions for 2x2 systems of hydrodynamic type in 
2 + 1 dimensions as obtained in |11| : 



uu 



tvir 



ipl + <ip u <p w - 2ip 

3 

V 



-<Pwip w + ipu^w - 2<fwVu + 2rj 

w^Pv 

V 



ipuw = , (36) 



r] w {<Pw ~ ipu) 

1 

7? 

1 

V 

vl. 

V 

we assume rj ^ 0: this is equivalent to the requirement that the dispersion relation of the system 
© defines an irreducible conic (indeed, the condition det(\A + B + riC) = is equivalent to A = 
ip+ip^+rjfi 2 ). We have verified that the system (f36l) is in involution, and all dispersionless limits 
appearing in Sect. 2.1 indeed satisfy these integrability conditions. Eqs. (|36|) are straightforward 
to solve. First of all, the equations for rj imply that, up to translations and rescalings, rj = 1, 
rj = u or rj = e w h(u). We will consider all three possibilities case-by-case below. Notice that ip 
and ip are defined up to additive constants which can always be set equal to zero via the Galilean 
transformations of the initial equation ([3]). Moreover, the system (|36p is form-invariant under 
transformations of the form 

Cp = ^ — sip + s 2 rj, ij) = ip — 2srj, rj = rj, u = it, w = w + su, (37) 

s=const, which correspond to the following transformations preserving the structure of equations 

x = x — sy, y = y, u = u, w = w + su. 
All our classification results are formulated modulo this equivalence. 



Vuu — 
Vuw — 

1]ww — 
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Case 1: rj = 1. Then the remaining equations imply ip = aw + f(u), ip = (3w + g(u), where / 
and g satisfy the linear ODEs 

f' = a(f'-P), g" = 2ag'-f3f'-(3 2 . 

The subcase a = leads to polynomial solutions of the form 

ip = Tit, <p = Pw- -(3(P + 7 )n 2 + 5u. (38) 
Up to equivalence transformations, the case leads to exponential solutions, 

if, = aw + 7 e an , <p = 5e 2au ; (39) 
here a, /3, 7, <5 are arbitrary constants. 

Case 2: r/ = u. Then the remaining equations imply ip = aw + f(u), (p = f3w + g(u), where / 
and g satisfy the linear ODEs 

uf" = a(f -(5)- 2/3, ug" = lag' - (if - (3 2 . 
The case a ^ {0,-1,-1/2} leads to power- like solutions of the form 

$ = aw + 7« a+1 , ip = 5u 2a+1 . (40) 
The subcase a = leads to logarithmic solutions, 

ip = -2/3ulnu- /Su, 99 = /3u> + P 2 uln 2 u + <5u. (41) 
The subcase a = —1 gives 

ip = — u; + 7lnn, ip = 5/u. (42) 
Finally, the subcase a = —1/2 gives 

■0 = — —w + 7-^1 99 = 51nu. (43) 

Case 3: rj = e w h(u). Then the remaining equations imply ip = e w f(u), ip = e w g(u) where /, g 
and h satisfy the nonlinear system of ODEs 

h" = f'-g, g"h = 2fg' -gf -g 2 , f"h = 2hg' - 2gti + ff - fg. 

Setting g = p' , f = h! + p, we can rewrite this system as a pair of third order ODEs 

hp'" = 2h'p" - p'h" + 2pp" - 2p' 2 , hti" = h'h" - 2h'p' + hp" + ph", 

which, up to a change of sign p — > —p, identically coincides with a system arising in the classifi- 
cation of integrable conservative hydrodynamic chains (subcase I\ of Sect. 3.1 in |13j). Setting 
p = h', the second equation will be satisfied identically, while the first one implies a fourth order 
ODE for h, h""h + 3(h") 2 - Ah'h'" = 0, whose general solution is an elliptic sigma-function: 
h = a(u), here (ln<r)" = —p, (p') 2 = 4p 3 — c (notice that <?2 =0, g% = c). Thus, as a particular 
case we have 

h = a(u), f = 2a'(u), g = a"(u). 
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Another subclass of solutions can be obtained by setting p = ch which implies 

h"'h - ti'ti = 2c(ti'h - h' 2 ) 

with the general solution 

h = ae (c+7)u + Pe (c ^ )u ; 

here a, f3, 7 are arbitrary constants. Although the structure of the general solution is quite 
complicated, one can show that Case 3 cannot arise as a dispersionless limit of an integrable 
third order soliton equation. 

3.2 Symmetric dispersionless limits 

In this section we consider first order equations of the form Q , 

U t = LpU x + t/jUy + TjWy + TV X , 
Vy — " X) 

where the coefficients ip, tp, rj, r are functions of u, v, w. We assume that the dispersion relation 
of this system defines an irreducible cubic, which is equivalent to the requirement 77 7^ and 
r 7^ (indeed, the dispersion relation has the form Xfi = r + <~p[i + tpfi 2 + t?/U 3 ). In this case 
the integrability conditions reduce to a system of first order partial differential equations for the 
coefficients ip, ip, rj and r which can be obtained from the requirement of the vanishing of the 
Haantjes tensor [12] as outlined in Sect. 3. The details are as follows: first we rewrite Eq. (jl]) 
in matrix form, 

Au t + Bu x + Cu y = 0, 
where u is a three-component column vector u = (n, v, w) 1 , and A, B, C are 3x3 matrices, 

/ -1 \ / ip t \ / v 77 \ 

A=\ , B= 1, C = -1 . 

\ooo/ V-100/ \oio/ 

The necessary conditions for integrability can be obtained from the requirement of the vanishing 
of the Haantjes tensor of the following family of matrices, 

{aA + (3B + 7 C)" 1 (aA + (3B + 7C). 

In fact, it is sufficient to require the vanishing of the Haantjes tensor for a two-parameter family 
(aA + B)~ l (aA + C). This condition turns out to be very restrictive, and leads to the following 
constraints for the coefficients (p, ip, r\ and r: 

T~u — ^Pv 1 

T~v — tyui 
77 

T w = 0, 

i>v = fw = 0, 



Vv = 0, 
- 1 

T 

Ttp w = ri<p v . 
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The integration of this system is straightforward. First of all, one can set ij) = f u , 77 = f w and 
'■P = 9ui T = 9v where / = f(u,w) and g = g(u,v). The separation of variables leads to the 
relations 

f w = a(w)k(u), g v = b(v)k(u), 

fuu = (3a(w)k(u), g uu = ab(v)k(u), 

where the functions a(w), b{v) and k(u) satisfy the ODEs a' = aa, b' = (3b and k" = a(3k; 
here a and (3 are arbitrary constants. Up to elementary translations, rescalings and Galilean 
transformations, this leads to the following subcases: 

Case 1. a = f3 = 0. This leads to equations of the form 

u t = v{uv) x + (J,(uw) y , 

where fx, v are arbitrary constants. These correspond to the Veselov-Novikov cases from Sect. 
2.2. 

Case 2. (3 = 0. This leads to equations of the form 

ut = v{uv + cra 3 /6) x + fi(e aw u) y , 

and 

u t = v(v + au 2 /2) x +n(e aw ) y , 
here fj,, u, a are arbitrary constants. 

Case 3. (3 7^ 0. This leads to equations of the form 

u t = u(e^k(u)) x + n(e aw k(u)) v , 
where v, jx, a, (3 are arbitrary constants, and k" = a(3k. 

4 Classification of integrable 3rd order dispersive equations 

Given an integrable dispersionless limit, one has to reconstruct dispersive terms. This can be 
done by requiring that all hydrodynamic reductions of the dipersionless system are inherited by 
its dispersive counterpart. We will illustrate this procedure using the KP equation, 

u t = uu x + w y + e 2 u xxx , w x = u y . 

Its dispersionless limit, the dKP equation, 

U t = UU X + Wy, W X = Uy, 

possesses one-phase solutions of the form u = R, w = w(R) where the phase R(x,y,t) satisfies 
a pair of Hopf-type equations 

R y = uR x , R t = (/x 2 + R)R X ; (44) 

here jx{R) is an arbitrary function, and w' = u. Equivalently, one can say that Eqs. (|44|) 
constitute a one-component hydrodynamic reduction of the dKP equation. Although the dKP 
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equation is known to possess infinitely many iV-component reductions for arbitrary N \15\ \TU[ 
\T7\ I18j. one-component reductions will be sufficient for our purposes. The main observation of 
[14] is that all one-component reductions (|44p can be deformed into reductions of the full KP 
equation by adding appropriate dispersive terms which are polynomial in the x-derivatives of R. 
Explicitly, one has the following formulae for the deformed one-phase solutions, 



u 



R, w = w(R) + e 2 L'R XX + hip." - {p'f)Rl \ + 0(e 4 ), (45) 



notice that one can always assume that u remains undeformed modulo the Miura group [7| . The 
deformed equations (|44|) take the form 



R y —\iR x 

+ e 2 Ur xx + h^" - {n'?)Rl\ + 0(6 4 ), 

Rt =(^ 2 + R)R X 

+ e 2 ((2/V + l)R xx + (////' - M (/x') 3 + (p') 2 mRl) x + 0(e 4 ). 



(46) 



In other words, the KP equation can be 'decoupled' into a pair of (1 + l)-dimensional equations 
(I46|) in infinitely many ways, indeed, n(R) is an arbitrary function. The series in (|45p and (|46j) 
contain only even powers of e, and do not terminate in general. 

Conversely, the requirement of the inheritance of all one-component reductions allows one 
to reconstruct dispersive terms: given the dKP equation, let us look for a third order dispersive 
extension in the form 

u t = uu x + Wy + e(...) + e 2 (...), w x = u y , (47) 

where the terms at e and e 2 are homogeneous differential polynomials in the x- and y-derivatives 
of u and w of the order two and three, respectively, whose coefficients are allowed to be arbitrary 
functions of u and w. We require that all one-component reductions (|44p can be deformed 
accordingly, so that we have the following analogues of Eqs. ([431) and (|46p. 

U = R, w = w(R) + e(...) + e 2 (...)+0(e 3 ), (48) 

and 

R y = nR x + e(...) + e 2 (...) + 0(e 3 ), R t = (/x 2 + R)R X + e(...) + e 2 (...) + 0(e 3 ), (49) 



respectively. In Eqs. (|48p and (|49p . dots denote terms which are polynomial in the derivatives 
of R. Substituting Eqs. (|4*8l) into (|47p . and using ([4^1) along with the consistency conditions 
= one arrives at a complicated set of relations allowing one to uniquely reconstruct 
dispersive terms in (|47[) : not surprisingly, we obtain that all terms at e vanish, while the terms 
at e 2 result in the familiar KP equation. Moreover, one only needs to perform calculations up 
to the order e 4 to arrive at this result! It is important to emphasize that the above procedure 
is required to work for arbitrary [i: whenever one obtains a differential polynomial in fi which 
has to vanish due to the consistency conditions, all its coefficients have to be set equal to zero 
independently. Another observation is that the reconstruction procedure does not necessarily 
lead to a unique dispersive extension as in the dKP case: one and the same dispersionless system 
may possess essentially non-equivalent dispersive extensions. In most of the cases one can get the 
necessary classification results working with one-component reductions only. There is however 
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one particular situation where one-component reductions are not sufficient. This is explained in 
the remark below. 

Remark 1. Let us consider the dKP equation, 

Ut = uu x + w y , w x = u y ; 

its one-component reductions (|44j) can be shown to satisfy a pair of additional first order con- 
straints, 

U 2 y - U X Wy = 0, (W t - UUy)U X - UyWy = 0. 

Conversely, any solution satisfying these constraints comes from one-component reductions. 
Similarly, one can show that two-component reductions of dKP are characterised by a pair of 
second order differential constraints, etc. Let us introduce an extension of dKP in the form 

u t = uu x + Wy + e(uy - u x w y ), w x = Uy\ 

by construction, it inherits all undeformed one-component reductions: the e-term vanishes on 
one-component reductions identically. This extension is, however, not integrable: one can show 
that it is not consistent with the requirement of the inheritance of TV-component reductions 
for N > 2. Thus, in what follows we eliminate deformations which inherit undeformed one- 
component reductions. 

In general, we proceed as follows. For definiteness, we will outline the algorithm for integrable 
dispersionless equations of the form ©, 

U t = ifU X + ifjUy + 7]Wy, W X = Uy. 

Its one-component reductions are of the form u = R, w = w(R) where R(x, y, t) satisfies a pair 
of Hopf-type equations 

R y = fiR x , R t = (p + ipn + r]n 2 )R x ; 

here n(R) is an arbitrary function, and w' = \x. We seek a third order dispersive deformation of 
Eq. (|3l) in the form 

u t = (pu x + ipu y + r]w y + e(...) + e 2 (...), w x = u y , 
and postulate that one-phase solutions can be deformed accordingly, 

u = R, w = w(R) + e(...) + e 2 (...) + 0(e 3 ), 

where 

R y = fiR x + e(...) + e 2 (...) + 0(e 3 ), R t = (if + # + W 2 )Rx + e(-) + e 2 (...) + 0(e 3 ). 

Proceeding as outlined above we reconstruct possible dispersive terms. In fact, one can start 
with arbitrary </>, ip, r/: our procedure will eventually recover the constraints obtained in Sect. 
3. However, using the classification results of Sect. 3 from the very beginning considerably 
simplifies the calculations. 

Remark 2. We point out that the formulae for dispersive deformations contain the expression 
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in the denominator. Since /i is assumed to be arbitrary, this expression is nonzero unless (p, tp, r] 
satisfy the relations 

Vw = 0, tp w + r] u = 0, (p w +ip u = 0, if u = 0. (50) 

These relations characterize the so-called totally linearly degenerate systems, which are known 
to be quite special from the point of view of the global existence of classical solutions: it was 
conjectured in [22] that smooth initial data for totally linearly degenerate systems do not break 
down in finite time. Modulo the integrability conditions (|36|) . the relations (|5U|) lead to equations 
of the form 

ut = a(wu x — uw x ) + (3(wu y — uw y ) + jw y , w x = u y , 

which have been discussed before in the context of the so-called 'universal hierarchy' [26J. For 
totally linearly degenerate systems (in particular, for linear systems), the procedure based on 
deformations of hydrodynamic reductions does not work, as the following simple example shows. 
Let us consider the KP equation, 

u t = auu x + w y + e 2 u xxx , w x = u y , 

where we introduced a parameter a: for a = the equation becomes linear. Looking for 
deformed one-phase solutions in the form 

u = R, w = w{R) +e 2 (...) +0(e 4 ), 

where 

Ry = nR x + e 2 (...) + 0(e 4 ), R t = (i? + aR)R x + e 2 (...) + 0(e 4 ), 

one can obtain the relation ab{R) — y! = where b(R) is the coefficient at R xxx in the e 2 -term 
in the expansion of R y . For a = one cannot solve for b(R), and obtains a relation fi' = 0. 
Thus, the linear equation ut = w y + e 2 u xxx does not inherit generic hydrodynamic reductions of 
its dispersionless limit. Another example of this kind is provided by the potential KP equation, 

u± = w y + |n 2 + e 2 u xxx . (51) 

One can show that this equation does not inherit hydrodynamic reductions of its dispersionless 
limit. However, some particular reductions can be inherited, for instance, those with /i=const. 

Thus, we exclude totally linearly degenerate systems from the further considerations: dis- 
persive deformations of such systems do not inherit hydrodynamic reductions, and require a 
different approach. 

4.1 Non-symmetric dispersive equations 

In this Section we summarize the classification results for integrable non-symmetric third order 
equations (H|), 

u t = tpu x + tpUy + f]w y + e(...) + e 2 (...), w x = u y , 

which are obtained by adding dispersive terms to integrable dispersionless candidates from Sect. 
3.1. Thus, we follow the classification of Sect. 3.1. 
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Case 1: We have verified that the exponential solutions (|39|) do not survive, so that all non- 
trivial examples come from the polynomial case (|38|) . 



(52) 



rj = 1, ip = ju, (p = /3w - ^/3(3 + y)u 2 + 5u. 

We point out that the corresponding dispersionless system possesses the Lax pair 
S y = 8uS x + r(S x ), 

St = \ 0w + + 7 )u 2 ) S x + duS x r'{S x ) + z^), 
where 

A 2/3+7 

r(^) = - 7 ^S. + S, /3 , z' = r' 2 . 

Lax pairs of this kind, consisting of two compatible Hamilton-Jacobi type equations, were first 
introduced by Zakharov in [35]. A detailed analysis of dispersive deformations leads to the two 
branches: 7 = 0, which corresponds to the (2 + l)-dimensional Gardner equation (Example 3 
of Sect. 2.1), and the case 7 = —3/3. In the latter case one can set 5 = 0, which leads to the 
apparently new equation (j5J), 

u t = (0w + 3 2 u 2 )u x - 3f3uu y + w y + e 2 [B 3 (u) - j3B 2 (u)u x ], 

where B = (3uD x — D y . The dispersionless limit of this equation possesses the Lax pair 

S x S v = BuS 2 H — , , , 

H x 3' (53) 

S t = f3 3 u 3 S 3 x -S3 + 3wS x , 
which follows from (I52h when 7 = —3/3. Its dispersive extension is 

l(j X y = BU11) XX + TTolP, 



3e 2T \ (54) 
yyy 



ipt = l3 3 e 2 u 3 tp xxx - e 2 ip yyy + 3p 2 e 2 uu y ip xx + f3wtp. 



This is case ([5]) from the Introduction. 

Case 2: One can prove that none of the logarithmic cases (|4ip . (|42p and ()43p survive, so that 
all non-trivial examples come from the power case (140 p . 

77 = u, ip = aw + ju a+1 , <f = 5u 2oi+1 . 

Further analysis leads to the following branches. 
Subcase 2.1: a = 1. In this case 

i] = u, i[) = w + yv 2 , <f = 5u 3 . 

The corresponding dispersionless Lax pair is of the form 

S y = ua, 

1 , ^ 3 ( 55 ) 
&t = uwa + -0(7 + a Ju , 
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where the function a(S x ) solves the ODE aa" — 2a' 2 = 35 + 27a'. The further analysis gives 
either 7 = 6 = 0, which leads to the non-symmetric Veselov-Novikov cases (Examples 4 and 5 
of Sect. 2.1, in this case one can take a = 1/S X ), or 5 = ^7 2 , in which case one arrives at the 
apparently new dispersive equation ([6]), 

u t = ^7 2 n 3 u a . + (w + ju 2 )u y + uw y + e 2 [B 3 (u) - ^u x B 2 (u)], 

where B = \^uD x + D y . This corresponds to the choice a = 1/S X — ^S x in the dispersionless 
Lax pair (|55l) , which gives 



Q Q 7 Q 2 u 

<->:r<->2/ — 7j z 

St = + S 3 y + ^u 3 S x + wS y . 

The dispersive extension of this Lax pair is 

7 1 

tpxy = ~^U^ XX - ^Uljj, 

e 2 7 3 o 2 f2 7 2 7 2 ^ 7 

Ipt = ~2J-U IpXXX + e Ipyyy —UU y 1p xx + —It a tj) x + 1^ ~ ^UUytp . 



(56) 



(57) 



The transformation 7 — ► 3/3, y — > —y, w — > —w reduces this case to Eq. ^ from the Introduc- 
tion. 

Subcase 2.2: a = —2. In this case one obtains 7 = 0, while 5 can be an arbitrary constant. The 
corresponding dispersive extension takes the form ([7]), 



u t = —^u x — 2wu v + uw v 

u 6 y y 



e 2 fl 



for 5 = it reduces to the Harry Dym equation (Example 6 of Sect. 2.1). The dispersionless 
limit of this equation possesses the Lax pair 



_ S x + r 

q 9 S * + T ■ ± S l + rS x 

S t = -2w V - 3 ; 

u z 3 u 6 

here r = 35/4. Its dispersive extension is of the form Lt = [A, L] where 



L "^ + 73^ + 4^' 

A _ 4g2 n3 , ( 6e2u * , 2 ^ew \ 2 5 ( 35u x yW 



(58) 



(59) 



Case 3: One can show that none of the examples from this class possess third order dispersive 
extensions. 
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4.2 Symmetric dispersive equations 

A detailed analysis of dispersive extensions of the form ([2]) 

u t = ipu x + ipu y + r]Wy + tv x + e(...) + e 2 (...) 



does not give any new examples: everything reduces to the two cases of Sect. 2.2. Notice 
that both symmetric VN and mVN equations can be viewed as linear combinations of the two 
commuting non-symmetric counterparts thereof. 

5 Concluding remarks 

We have proposed a new approach to the classification of integrable equations in 2+1 dimensions 
based on the concept of hydrodynamic reductions and their dispersive deformations. It consists 
of the two steps: 

— Classification of dispersionless systems which may (potentially) arise as dispersionless limits of 
soliton equations. This can be efficiently achieved using the method of hydrodynamic reductions 
as outlined in [TOl: 



— Classification of possible dispersive deformations based on the requirement that hydrodynamic 
reductions of the dispersionless limit are inherited by the dispersive equation |14j . 
This procedure was applied to the classification of third order soliton equations with 'simplest' 
nonlocalities. Further research in this direction may include the following topics: 

(a) Classification of more general (in particular, higher order) soliton equations/systems with 
more complicated structure of nonlocal terms. Thus, one may allow 'nested' nonlocalities of the 
type w = D~ l D y u, v = D^ 1 D y F(u,w), etc. 

(b) Construction of dispersive deformations via an appropriate quantization of the corresponding 
dispersionless Lax pairs [55] . 

(c) Investigation of the structure of multi-soliton solutions of the new equations ([5]) - ([7|) in the 
spirit of [3 [6] . 
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